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Abstract 

Paramagnetic pair breaking is believed to be of increasing importance in 
many layered superconducting materials such as cuprates and organic com- 
pounds. Recently, strong evidence for a phase transition to the Fulde-Ferrell- 
Larkin-Ovchinnikov(FFLO) state has been obtained for the first time. We 
present a new theory of competing spin and orbital pair breaking in clean su- 
perconducting films or layers. As a general result, we find that the influence 
of orbital pair breaking on the paramagnetically limited phase boundary is 
rather strong, and its neglect seldom justified. This is particularly true for 
the FFLO state which can be destroyed by a very small orbital contribution. 
We discuss the situation in YBa2Cu3C>7 which has two coupled conducting 
Cu-0 layers per unit cell. As a consequence, an intrinsic orbital pair breaking 

component might exist even for applied field exactly parallel to the layers. 
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I. INTRODUCTION 



Most theoretical studies of paramagnetic pair breaking in superconductors followed the 
attitude of the classical papers by ClogstonEl and ChandrasekhaJl where only spin pair 
breaking was considered and the orbital component was assumed to be negligibly small. A 
notable exception is the dirty limit theory developed by Makil, FuldJ and others. In many 
experiments, on the other hand, both pair breaking components are present and the neglect of 
the orbital contribution is not really justified. Recently, ultra-thin films became available and 
several new classes of layered superconducting compounds have been discovered. For applied 
field parallel to the filmsil or conducting planes!, Pauli paramagnetism can be the dominating 
pair breaking effect, provided the conducting layers are sufficiently separated from each other 
or the thickness of the films is sufficiently small. In many of these compounds, including 
High-T c cuprates and organic superconductors, impurity scattering and spin-orbit coupling 
is small and orbital pair breaking is - for an applied field parallel to the planes - the most 
important second order effect, next to the spin effect, to be taken into account. 

Of particular interest is the Fulde-Ferrell-Larkin-Ovchinnikov(FFLO) state&i, which is a 
spatially inhomogeneous superconducting state, predicted to occur in clean superconductors 
with purely paramagnetic limiting. Recent critical field measurements! in the quasi-two- 
dimensional organic superconductor k — (BEDT-TTF) 2 Cu(NCS) 2 strongly suggest that a 
state of the FFLO type exists in this material; agreement between experiment^ and existing 
theories has been successfully checked^ both in view of the angle- dependence^ and the 
temperature dependence^ of the upper critical field (see alsolll). Apparently, this is the first 
time since the original predictions in 196401 that quantitative agreement between theory 
and experiment with regard to the FFLO phase boundary has been established. Strong 
paramagnetic effects can also be expected for the High-T c cuprate superconductors at low 
temperatures, when the conducting planes in adjacent unit cells are well separated from 
each other. A measurementi'0 at T = 1.6 K in YBa2Cus07 indicates rather clearly that 
the superconducting state is paramagnetically limited but, on the other hand, the observed 
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transition is too broad to allow a decision between the FFLO state and the homogeneous 
superconducting state. 

A measure of the relative strength of orbital and paramagnetic pair breaking is the ratio 
of the paramagnetic critical field H p divided by the orbital upper critical field H c2 of a type 
II superconductor. For a bulk superconductor in the clean limit this ratio can be written as 

H P /H c2 ~ & (1) 

Kp 

in terms of the Fermi wavelength kp and the coherence length £o of BCS theory. 

This relation implies that orbital pair breaking will always be the dominating mechanism 
in bulk superconductors, no matter how large the Ginzburg-Landau (GL) parameter k is; 
this holds at least in the framework of conventional BCS theory. In a thin superconducting 
layer of thickness d < £o , on the other hand, the orbital critical field H c (d) is increased by 
a factor of £o/d and the corresponding ratio is given by 

H p /H c (d) ~ -i. (2) 

Kp 

In comparison to Eq. (JJ) a small transverse dimension d <C £o of the film suppresses the 
orbital effect and enlarges the spin effect drastically. However, equation (0) also shows, 
that the critical thickness which separates the spin pair-breaking and orbital pair-breaking 
dominated regimes is still of the order of an atomic distance. Thus, the estimate (Q), which 
is confirmed by more quantitative calculations to be presented below, indicates that a nearly 
perfect two-dimensional situation is required in order to justify the neglect of orbital pair 
breaking contribution in clean superconductors. (The situation in dirty superconductors 
is much more favorable for the spin effect; the FFLO state, however, is suppressed by 
impurities). 

The simultaneous action of both types of pair breaking has already been studied for a 
particular situation, an infinitely thin superconducting film in a tilted magnetic field00@. 
In such a configuration, orbital pair breaking is entirely due to the perpendicular field 
component, while the component parallel to the film is exclusively responsible for the spin 



effect. The states near the upper critical field can be characterized by different Landau 
quantum numbers n = 0,1,2,..., depending on the tilt angleElE. The order parameter 
structure near H C 2 has also been investigated and shows interesting properties, such as 
several zeros with different vorticity per unit celfi For exactly plane-parallel external field 
the FFLO phase is recovered. However, in this limit no orbital pair breaking exists and only 
the spin effect survives, as a consequence of the vanishing thickness of the superconducting 
layer in this model. 

In this paper we investigate a superconducting film of finite thickness in a magnetic 
field parallel to the conducting plane. Thus, the usual model of purely paramagnetic pair 
breaking is generalized in a different way, taking into account the influence of a finite orbital 
pair breaking component on the FFLO state. The model is formulated in section ||, using 
the framework of the quasiclassical Eilenberger equations. We assume that the film thickness 
is smaller than the coherence length and use a cylindrical Fermi surface. This shape of the 
Fermi surface allows us to study the influence of orbital pair breaking without any additional 
complications, like scattering of quasiparticles at the film boundaries. Such boundary effects 
seem less important in the present context, but may, nevertheless, be present in many 
materials and should be taken into account in future work. In addition, the cylindrical 
shape of the Fermi surface, which corresponds to a truly two-dimensional situation, allows 
us to extend our investigations to superconducting layers of atomic dimension. The main 
results, obtained by solving numerically the relevant phase boundary and stability equations, 
and a discussion of possible orbital pair breaking contributions in the plane-parallel field 
configuration of YBa2Cu307, are reported in section pTj. Finally, the results are summarized 



in section IV 



II. BASIC EQUATIONS 

We first calculate the highest field where a superconducting solution of the quasiclassical 
equations, for small order parameter A, exists in a thin film. This field may correspond 
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to a second order phase transition or to the supercooling limit of the normal conducting 
state; to make a decision between these two possibilities the free energy of the competing 
homogeneous superconducting state will be calculated in a second step. 



A. Stability limit of normal conducting state 

Let the film be parallel to the xy plane with a finite extension from — d/2 to +d/2 in 
the z— direction. The applied magnetic field B is assumed to be parallel to the plane of the 
film and parallel to the y— direction, B = Be y . The transport equations, linearized in A, 
are given by 

2(w, - ifiB) + M F {k) (V r - i{2e/hc)A)] f(f, k, u a ) = 2A(r, k), (3) 

2(ui - i/iB) - M F (k) (V r + i(2e/hc)£)] f + {f, k, u s ) = 2A*(f, k), (4) 

Here, the Zeeman term [iB occurs in the combination uj s = ooi—ifiB, where u>i = (2/+l)7r/csT 
are Matsubara frequencies, /i ~ h\e\/(2mc) is the magnetic moment of the electron and B 
is the magnitude of the induction. The self-consistency equation for the gap is given by 

Nd 1 t\ n d 
27rk B T^- + ln- \A(f,k)=irk B T J £ i d 2 k'V(k,k') \f(r,k> > s ) + f(r,k' >,*) , (5) 

where Nd is the cutoff index for the Matsubara sums. The Fermi velocity is given by vp{k) = 
vp(e x cosip + e^sin^) = vpk\ the integral in Eq. (|5p over the cylindrical Fermi surface is 
simply a one-dimensional integral over the angle variable (p. We allow for a separable gap 
anisotropy, A(f,k) = A(r)y(k), V(k,k') = y(k)"f(k'), which will be specialized later to s- 
wave and d-wave superconductivity. We use the following gauge for the vector potential: 
A x = Bz, and A y = A z = 0. 

The standard method to solve the linearized transport equations uses a complete set 
of eigenf unctions of the operator kd r to construct the inverse of the differential operators 
on the l.h.s. of Eqs. ©, (H). Here, d r is an abbreviation for the gauge- invariant derivative, 
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d r = V r — i(2e/hc)A. For a cylindrical Fermi surface, kd r contains no derivative with respect 
to z and the Green's functions depend on z in a purely local way (z playing the role of a 
parameter). This allows a straightforward generalization of the standard method to the 
present problem. 

Let us start from the well-known bulk solution of eqs. If the eigenf unctions and 

eigenvalues of the operator kd T for an infinite sample are denoted by fy and ikp respectively, 

k ^r fijf) = ikp fajf) , (6) 
then the solution of Eq. @ is given by 

f(r,k,u s ) = dn z——A(ri,k), (7) 

J J [Ztt) 6 Us + iv F kp/2 

where the spatial integration extends over all space. In the chosen gauge the solutions of 
Eq. (^|) are given by 

with the abbreviation K\\ = ^13. Using the completeness of the set of eigenfunctions @ 
the Green's functions may immediately be written in the form of Eq. (0). 

To transform relation (|7|) to a finite volume, it is, in our case, only necessary to restrict 
the spatial integration in Eq. (0) to the film volume, i.e. to perform the integration over z\ 
from — d/2 to +d/2. This simple method works only for a cylindrical Fermi surface, where 
the momentum of the quasiparticles is always parallel to the film boundaries. Otherwise, 
quasiparticle scattering at the film boundaries leads, for small d < £q, to a modification of 
the integral kernel which has to be calculated by solving Eq. (||) in a finite volume, with 
appropriate boundary conditions. 

To proceed, the denominator of the integrand in Eq. (|7|) is shifted into an argument of 
an exponential function by means of the identity 
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dte~ rt . (9) 



Now, if the eigenfunctions (|8|) are inserted, the Green's function (|7|) is represented as an 
integral with regard to the variables p, f, and t over an exponential function. Two of these 
integrations can be performed analytically and the Green's function takes the form 

POO ~ /• A A 

f(r,k,cu s ) = dte- tuJ °e—z tVFK " zk *A(r--v F k,k). (10) 
Jo 2 

For d < £ the order parameter may be considered as ^-independent and the Greens's 
function f(f, k, tu s ) may be replaced by its value f(f, k,cu s ) which depends only on x and y 
and denotes the average of f(f, k, ou s ), with respect to z, from — d/2 to +d/2: 

f°° 1 /l ~ \ t - * 

f(r,k,u a )= dte~ tuJ °- ^-sin -t^K,^ A(f - -v F fc, fc). (11) 

The rest of the calculation is a straightforward generalization of methods developed in 
previous worksite. It is convenient to perform the following shift in the argument of the 
space-dependent part of the gap: 

A(f - —vpk) = e-5 to ^A(f). (12) 

Inserting the Green's function solutions in the self-consistency equation for the gap (^) 
and performing the Matsubara sum yields the linearized gap equation 
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-ln('^ A(f)=nk B T f° dt— ^— - t —l{k'f 
\TJ y ' Jo sinh(vrTt) Jo 2vr ,v ; 

1 .*,\ ( „. 1 
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sin ^~tvFK\\dk'^J cos {^iBt — —tvpk'V 



\tv F K\\dk' x 
A(f) (13) 



which has to be solved in order to find the magnetic field where the normal-conducting state 
breaks down. The operator V used in Eq. (|I3"D acts in the x, |/-plane. Eq. (p~3f ) differs from 
previous results^ by a <i-dependent factor which reduces to 1 in the limit d — > 0. For d-wave 
superconductivity the finite thickness of the film breaks rotational invariance; if $ is the 
angle between the magnetic field and the y— axis of the crystal, the following replacement 
has to be performed in the integrand of Eq. fli"3]): 



dk' x =^ d (k' x cos<& — ^sin$) , (14) 



in order to take the angle dependence of the external field into account. 

To proceed further, the gap A is assumed to be proportional to plane wave states e^ r . 
Solving the linearized gap equation (|T3|), with V replaced by iq, for different wave numbers 
q, one obtains a function B(q). The field we are looking for, where the normal conducting 
state breaks down - and the corresponding wave number - is given by the highest B(q). 

Eq. (|l^) comprises two limiting cases where the behavior of the solutions is known, the 
purely paramagnetic limit (for d = 0), and the purely orbital limit (for /x = 0). For d = 
one obtains the standard FFLO resultS,!. For /i = Eq. (fl3|) may be solved analytically 
near T c . In this limit one obtains for the parallel critical field, where the normal- conducting 
solution breaks down 
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Eq. (|15D is in agreement with the second order transition line obtained for d < £(T) in 
the framework of GL theory. At low temperatures microscopic calculations have, to our 
knowledge, only be performed for Fermi surfaces of spherical0^0 or ellipsoidal^ shape. In 
these works a d~ 3 ^ 2 behavior for the critical field has been obtained in the limit of very small 
d. In contrast, the present theory, using a cylindrical Fermi surface, leads to an approximate 
d^ 1 behavior of the critical field in the whole temperature range. 



B. Phase boundaries of homogeneous states 

For an infinitely thin film, the free energy of the FFLO state is only slightly lower than 
the free energy of the homogeneous (paramagnetically limited) superconducting state. The 
presence of an orbital pair breaking component, realized by a vector potential, in our film 
of thickness d may change the free energy balance in a decisive way. To clarify this point, 
we calculate the free energies of the homogeneous superconducting and normal-conducting 
states and compare the resulting phase boundaries and stability limits with the FFLO tran- 
sition line studied in the last subsection. Considering films or layers of finite thickness, we 
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refer to states which do not depend on the coordinates x, y within the plane as 'homoge- 
neous states'; these states may, nevertheless, depend on z as a consequence of a (residual) 
screening property of the thin film. 

The spatially constant, purely paramagnetically limited superconducting state has been 
studied first by Sarma@; this case corresponds to d = in the present model. For T/T c < 
0.56 he found a first order phase boundary, which lies below the FFLO transition line. This 
line is determined by inserting the solutions of the nonlinear gap equation 

N D 1 rp N ° ( 1 \ 

27rkB Tj2- + ^- = 7rk B TY: , +cc , (16) 



into the free energy difference 
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F s - F n = -N(E F )nk B TY: ± / J + cc . (17) 



~ \ VI A I 2 + ^ 2 

For T < 0.56 T c the function of B has two branchesSH, as shown in Fig. [L| Our 

calculation of the second variation of the free energy shows (see the curves for T/T c = 0.1 in 
Fig. [I]) that the homogeneous superconducting state may be superheated up to the highest 
field, where the two branches cross. The lowest field where the lower branch exists defines, 
on the other hand, the supercooling limit of the normal-conducting state. Thus, the region 
of the lower branch corresponds, as expected, exactly to the metastable region of the first 
order transition. In this way, three transition lines, the phase transition line where the free 
energies coincide, the superheating line, and the supercooling line, are determined by solving 
Eqs. dUD, (0); for finite d the same method is used to determine the metastable region. 

For the present circular Fermi surface, the determination of the 'homogeneous states' 
in a film of finite thickness d is still a local problem despite the nontrivial z— dependence 
appearing in the transport equations. The assumption of a gap which depends weakly on 
the z— coordinate leads, in analogy to the reasoning of the last subsection, to the following 
relation between the averaged Green's functions /, g and the gap A: 

/=_i_arctaa(^))A^ (18) 
B{d,tp) lo s 



where 



B(d,<p) 



VpKndcOSlfi 



(19) 



47rk B T c 

Note, that orbital pair-breaking leads to a dependence of the Green's functions on the 
quasiparticle wave number k. Using Eq. flI%D the self-consistency relation for the gap takes 
the form 
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+ cc. 



(20) 



and the free energy difference is given by 



F s - F n = —N(E F )-nk B T ^ [** ^ 

Tr n JO Z7T 



V 



No r2n 

E 

(=0 

where the factor Ai(d, <p) is defined by 

1 1 



(Md, <p) - yJiAmtf + MdM' 



+ cc. 



(21) 



,B(d, ip). 
arctanf ) . 



(22) 



A t (d,<p) B(d,<p) 

The factor P(<p) is 1 for s-wave and l + cos(4<^) for d-wave superconductivity. Eqs. (p0|) , (pTf) 
are essentially of the same (local) form as Eqs. (|T6f) , (|T7|) ; the main difference is the depen- 
dence of the Green's functions on the direction (p of the quasiparticle momentum, which 
leads to the ^—integrals in Eqs. (pOD, (pTD- For d — > the previous results are recovered as 
Ai(d, (p) approaches uj s in this limit. Solving Eq. (p0[) and calculating the free energy ( |2T|) 
the influence of a finite orbital pair breaking contribution, due to a nonzero film thickness 
d, on the three transition lines can be studied. 



III. RESULTS AND DISCUSSION 

In this section we discuss the following four transition lines, defined in more detail in the 
last section: (1) the line B c where the free energies of the homogeneous superconducting and 
normal-conducting states coincide, (2) the superheating limit B s h of the homogeneous super- 
conducting state, (3) the supercooling limit (stability limit) of the normal-conducting state 
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against spatially homogeneous superconducting fluctuations, which is denoted by B sc , and 
(4) the stability limit of the normal-conducting state against spatially inhomogeneous super- 
conducting fluctuations, denoted by Bfflo- The detailed results reported here are mainly 
for isotropic gap (s-wave) superconductors; a few calculations have also been performed for 
d-wave superconductors in view of recent experiments on YBa 2 Cu30y. 

As a starting point, we show in Fig. § these four transition lines for d = 0, in the purely 
paramagnetic limit. For T > T tr i = 0.56 T c all four lines merge into a single second order 
transition line. Below the tricritical point T tr i only three different lines are visible in Fig. |2| 
since, interestingly, Bfflo and B s h exactly coincide (this coincidence occurs, however, only 
for a circular Fermi surface). Only Bfflo is physically significant for d = 0, the other lines 
are meaningless. At Bfflo a second order phase transition to the FFLO state takes place; 
the lower phase boundary of the FFLO state is not dealt with here; it has been studied by 
Burkhardt and Rainer0. 

The paramagnetic pair-breaking effect dominates for very small d while the orbital effect 
dominates for large d. Thus, it should be possible to define a critical thickness d* which 
roughly separates the two regimes. This crossover behavior is shown in Fig. ^ for the 
thermodynamic critical field B c . The value of d* is of the order of k F l in the region of low 
T, in agreement with the estimate of section |. Fig. ||| shows also the decreasing importance 
of paramagnetic pair breaking with increasing T. Generally, the additional orbital pair- 
breaking effect brought about by the finite thickness of the conducting layer, leads to a 
depression of all four fields shown in the reference figure (Fig. |2|). A detailed plot of the 
T— dependence of the fields B F flo, B c , B sc (the superheating field B sh lies above B F flo an d 
has been omitted for clarity) for d/k F l = 0.5, 1.0, 2.5 as shown in Fig. f| reveals, however, 
significant differences. The FFLO transition field is much stronger suppressed than the 
line B c where the free energies of the homogeneous states coincide. As a consequence, 
for d ^ 0.5 kp 1 (see top of Fig. |) the FFLO transition vanishes {Bfflo becomes the 
supercooling field of the normal state) and is replaced by a first order transition at B c to 
the homogeneous superconducting state. Note that a conducting layer of atomic thickness 
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(dimension of unit cell in the plane) yields enough orbital pair breaking to produce this 
suppression of B FF lq in favor of B c . This behavior is not unreasonable; spatially varying 
states are known to be much more sensitive to perturbations than homogeneous states 
(recall in this context Andersons theorem). With further increasing orbital pair breaking 
(see middle and bottom of Fig. [|) the lines B FFLO and B sc tend to merge and the metastable 
region shrinks; for d > 3 kp 1 orbital pair-breaking dominates. 

The behavior of all four fields as a function of d is shown in Fig. [| in the low temperature 
region (for T = 0.01 T c ), where paramagnetic effects are most pronounced. This figure gives 
an overview of the cross-over from the paramagnetically dominated regime at very small d 
to the orbitally dominated regime at large d. If the thickness where B FFLO and B c cross is 
denoted by di, then the FFLO state is only realized in the small range d < d\ = 0.5 k F , 
for d > d\ a first order transition to the homogeneous (mainly) paramagnetically limited 
state occurs. The FFLO line plays the role of a supercooling limit of the normal state until 
it falls (at d = 1.2k F ) below the line B sc , where the normal- conducting state is limited 
by spatially constant superconducting fluctuations. The wavenumber q of the FFLO state 
decreases with increasing d until it jumps to q = at the line B sc (at the crossing point two 
degenerate solutions exist for q). It should be pointed out that we continue to use here the 
term 'FFLO state', even if this term denotes, strictly speaking, a state without any orbital 
pair breaking contribution. With increasing d the lines B s h, B c , B sc approach each other 
and the transition becomes identical to the well-known second order transition of a thin film 
in a parallel field, which is entirely due to orbital pair breaking (the region of really large 
d where the difference between type I and type II superconductivity becomes important is 
clearly outside the range of validity of the present model). 

The results reported so far have been restricted to s-wave superconductors, with isotropic 
gap [using 7 (A;) = 1 in Eqs. (|1"3"D, (|20|) ] . A few calculations have also been performed for d- 
wave superconductors, using 7 (A:) = 2(A^ — ky). The results generally confirm the behavior 
found for s-wave superconductors. An interesting peculiarity of d-wave superconductors 
without any orbital pair-breaking is a steep rise of B FF lo with decreasing T below T/T c = 
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0.1 (see Fig. 4 of Manalo and Kleinc2l). This peak belongs to the <p q = portion of the 
critical fieldi and is much steeper than the corresponding part of the critical field curve for 
s-wave superconductors. Our calculations show that this peak can be effectively suppressed 
by a very small (d/kp 1 « 0.2) amount of orbital pair breaking. The absence of this peak 
in measurement M on k — (BEDT-TTF) 2 Cu(NCS) 2 may be an indication of a very small 
residual orbital pair-breaking contribution in this material; a more detailed study would be 
worth while. 

Let us investigate the consequences of the fact that the FFLO state can be suppressed in 
favor of the homogeneous superconducting state by a very small admixture of orbital pair- 
breaking (see Figure ||) . We have not mapped out a complete phase diagram like figure |5] for 
d-wave superconductors (where B FFLO becomes anisotropic as a consequence of the finite d) 
but instead performed a few calculations in order to get an overview of what happens. The 
results confirm qualitatively the main feature visible in figure |5|, namely a much stronger 
suppression of B FFLO , as compared to B c , by orbital pair-breaking. For YBa 2 Cu 3 07 we have 
two conducting Cu02 layers per unit cell with a distance of ~ 3.9 A, while the c-axis zero- 
temperature coherence length is estimatedil to be 2 — 4 A. The coupling between bilayers 
in adjacent unit cells may obviously be neglected, as a consequence of the large length 
c ~ 11.7 A of the unit cell in this direction. The coupling between the two layers in one unit 
cell, on the other hand, remains an open question, and the following two possibilities should 
be taken into consideration. 

The first possibility is, that the two superconducting layers decouple at low T, below 
some crossover temperature T*. As is well known, the orbital critical field of weakly coupled 
layers diverges^ below some crossover temperature T*, which means that paramagnetic pair- 
breaking is the only remaining mechanism to limit the superconducting state. This requires 
a two-dimensional, in-plane mechanism of superconductivity. The amount of orbital pair- 
breaking would be negligibly small in this case and the superconducting state below the 
critical field should be the FFLO state. 

The second possibility is that the superconducting state keeps its finite extension for 
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arbitrary T. This requires an inter-plane mechanism where the bilayer structure is essential 
for the superconducting pairing process. In this case, the bilayer may be approximately 
replaced by a single layer of finite thickness d ^ 2 — AA . Taking a value of Vp ~ 7 ■ 
10 7 cm/sec for the Fermi velocity in the a — b plane, as measured by Andreev reflections^, 
we estimate a value between 1 and 4 for our dimensionless thickness parameter d/kp 1 , which 
measures the amount of orbital pair-breaking. Thus in this case, if the bilayer structure can 
be approximated by a finite slab, the amount of intrinsic orbital pair-breaking in YBCO, 
brought about by the finite thickness of this slab, will be by far large enough to suppress the 
FFLO state. The second order FFLO transiton will be replaced by a first order transition to 
a homogeneous superconducting state; this transition is due to the combined action of both 
pair breaking mechanisms rather than a single one. For a typical value of d/kp 1 « 2, the 
critical field would be still of the same order of magnitude as the purely paramagnetically 
limited, (Pauli limiting) field at d = but with a strongly reduced metastability region (see 
Fig-D 

The transport measurement of Dzurak et al.i of the critical field of YP^CusCV at 1.6K 
led to a result of the order of the Pauli limiting field B p = B c (d = 0); the transition seems, 
however, too broad to distinguish between the FFLO and Pauli limiting fields. Further, 
more accurate experiments are required to settle this question, which concerns fundamental 
aspects of the superconducting state in High-T c cuprates. Observation of the FFLO state in 
the plane-parallel field configuration of YBa 2 Cu307 would be a strong argument in favor of 
an in-plane mechanism of superconductivity in this material. The relevance of this question 
has also been discussed by Yang and Sondhfi using the framework of the Lawrence- Doniach 
model, which is in a sense complementary to the present approach. 

The orbital pair-breaking effect due to the finite thickness of the conducting layers has 
previously been taken into account in a theory by Schneider and Schmidt@. This theory may 
be used successfully to fit the experimental data0 near T c , but neglects all paramagnetic 
effects. The purely paramagnetic limit, on the other hand, has been studied by Maki and 
Won@ and by Yang and Sondhi0. The present results show that both effects should be 
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taken into account for a detailed description of the transition. The orbital effect cannot 
be neglected even if the orbital critical field is several times higher than the paramagnetic 
limiting field. 

Finally, we note, that single atomic (molecular) layers are responsible for the supercon- 
ducting state in the organic compound k — (BEDT-TTF) 2 Cu(NCS) 2 where phase boundaries 
compatible with a d-wave version of the FFLO state have recently been observed. Thus, 
a considerable influence of orbital pair-breaking, which would suppress the FFLO state in 
favor of the homogeneous superconducting state, does not exist in this material. Measure- 
ments of the detailed temperature-dependence of Bfflo close to T = could reveal even 
small admixtures of an orbital pair-breaking component. 

IV. CONCLUSION 

We developed a theory of competing paramagnetic and orbital pair-breaking effects in 
clean superconducting films and layers. The destructive influence of orbital pair-breaking on 
the FFLO state turned out to be stronger than commonly expected. It is necessary to have 
single-atomic layers in exactly plane-parallel fields in order to be able to neglect completely 
the orbital component and observe the "pure" FFLO state. We calculated only the upper 
phase boundary of the "mixed" FFLO state (which occurs in films of finite thickness as a 
consequence of the combined action of both pair breaking mechanisms). The equilibrium 
structure as well as the lower transition line to the homogeneous superconducting state, 
have not yet been calculated for this new inhomogeneous state. Another challenging open 
question is a microscopic treatment of Josephson coupling between layers of finite thickness. 
Our results show that a careful measurement of the plane-parallel critical field in YBa 2 Cu 3 07 
could give important information on the question whether an in-plane or an inter-plane 
mechanism is responsible for superconductivity in High-T c superconductors. 
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FIGURES 

FIG. 1. Stable and unstable branch of A as a function of B for t = T/T c = 0.1 and second 
derivative of free energy Fa* a for t = 0.1, 0.3, 0.6. 

FIG. 2. Critical fields Bfflo-, B c , and B sc as a function of T for d = 0. The superheating field 
B s h agrees with Bfflo for the present circular Fermi surface. 

FIG. 3. Critical field B c as a function of the dimensionless thickness parameter kp 1 jd for 
t = T/T c = 0.5 and t = 0.9. The purely orbital pairbreaking curve, which is proportional to d^ 1 
and the purely paramagnetic curve, which is independent of d are shown as dotted lines (t = 0.5). 

FIG. 4. Critical fields B FFLO , B c , B sc as a function of T for d/kp 1 — 0.5 (top), d/kp 1 — 1.0 
(middle), und d/kp 1 = 2.5 (bottom). 



FIG. 5. Critical fields Bfflo, B c , B sc , B s h as a function of d for T = 0.01 T c . 
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